Gaussian potentials facilitate access to quantum Hall states in rotating Bose gases 
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Through exact numerical diagonalization for small numbers of atoms, we show that it is possible 
to access quantum Hall states in harmonically confined Bose gases at rotation frequencies well 
below the centrifugal limit by applying a repulsive Gaussian potential at the trap center. The main 
idea is to reduce or eliminate the effective trapping frequency in regions where the particle density 
is appreciable. The critical rotation frequency required to obtain the bosonic Laughlin state can 
be fixed at an experimentally accessible value by choosing an applied Gaussian whose amplitude 
increases linearly with the number of atoms while its width increases as the square root. 
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For almost 30 years since its initial observation in a 
two-dimensional electron gas [1], the quantum Hall (QH) 
effect has continued to be the subject of intense research. 
Much recent theoretical work indicates that similar in- 
compressible states should also exist in ultracold trapped 
atomic gases [2, 3]; the most intriguing possibility would 
be neutral bosons with tunable two-particle interactions. 
Many of the proposals exploit the formal equivalence of 
the Hamiltonian for a two-dimensional interacting elec- 
tron gas in the presence of an external transverse mag- 
netic field with that describing a neutral interacting gas 
subject to rotation [4, 5]. Other theoretical approaches 
make use of internal degrees of freedom to create artificial 
magnetic fields [6-9]. 

In addition to their inherent interest, fractional QH 
states may also be used to perform intrinsically fault- 
tolerant quantum computing. The topological properties 
of quasiparticle excitations with non-Abelian statistics 
could be used to generate a universal set of multi-qubit 
gates that are protected against the deleterious effects 
of decoherence [10] by braiding them around each other 
in specific patterns [11]. The QH state with filling fac- 
tor v = 5/2 in electronic systems is widely believed to 
be described by a Moore-Read wavefunction which pos- 
sesses the required properties, though recent work has 
cast some doubt on this [12, 13]. The equivalent state 
in rotating Bose gases (RBGs) is expected at = 1 [5]. 
A distinct advantage of ultracold atomic gases for topo- 
logical quantum computing is that quasiparticle creation 
and manipulation could be as straightforward as shining 
blue-detuned lasers into the gas and moving them around 
each other [14]. 

Unfortunately, the QH effect has not yet been observed 
in RBGs. Incompressible states result when the number 
of quantized vortices (the analog of quantized flux in 
electronic systems) is comparable to the number of par- 
ticles N, so that the filling factor v = N/N^ ~ 1. This 
strongly correlated regime generally requires either very 
low particle number, or very high rotation frequencies 
approaching the centrifugal limit of the confining har- 
monic trap [15]. The first issue can be addressed by 



adding a deep one-dimensional optical lattice oriented 
along the axis of rotation, making an array of discon- 
nected pancake-shaped wells [16], each containing a small 
number of atoms (conceivably on the order of lO'^), of 
which all but one or two can be eliminated [17]. The 
second issue can be partly resolved by applying an ap- 
proximately quartic potential to the harmonic trap (in 
practice by adding a repulsive Gaussian potential) [18], 
though at high rotation frequencies this appears to favor 
weakly correlated 'giant vortex' states [19]. 

Using exact calculations for small numbers of parti- 
cles, we demonstrate that the application of a repulsive 
Gaussian potential with the appropriate parameters can 
yield bosonic fractional quantum Hall states in RBGs at 
experimentally accessible rotation frequencies. A simple 
model agrees closely with the numerical results. To favor 
the Laughlin state, the combined potentials (harmonic 
trap, centrifugal term, and external Gaussian) should be 
flat over the spatial region where the particle density is 
appreciable. This requires a Gaussian whose amplitude 
increases linearly with the number of bosons N (with a 
small coefficient) while the width scales as \/N. 



The Bose gas at zero temperature is harmonically 
trapped in an axisymmetric potential with radial and 
axial trapping frequencies lo and ujz, respectively, and 
is rotated around the z-axis at a rate fl. Imposing 
tight confinement along the z-axis (tUz ^ w), all motion 
along this direction is frozen out and the cloud is quasi- 
two-dimensional, with each particle's position denoted 
Ti = {pi, 4>i)- An external Gaussian potential of the form 



76" 



is included to represent a repulsive blue- 



detuned laser that is being shined upon the center of the 
rotating Bose gas, with tunable amplitude 7 and width 
a. The bosons interact via the usual delta- function pseu- 
dopotential with a two-dimensional coupling constant 
g — V87r?iw^^a/£z, where variables t — yJh/Muj and 
— \/1T'/Mu>z are the characteristic oscillator lengths 
along the radial and axial directions, respectively, and a 
is the three-dimensional scattering length. In the rotat- 
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FIG. 1: The yrast lines for the six-particle RBG for different 
external potential intensities 7. Parameter values are a = 2.0 
and g = 0.1. 

ing frame the Hamiltonian is thus 



N 



H 



£ \^ - Mf2 X r,)' + ^u^lzf (1) 



i=l 

M 



where A'' is the number of bosons and M their mass. 
The energy spectrum of the RBG can then be calculated 
by exact diagonalization of the Hamiltonian H, for var- 
ious values of the parameters Q = Q/oj, g = g/Utji'^, 
7 = ^//h,uj and (J = a 1 1. For more details about the 
diagonalization process, see Ref. [20] . 

A Bose condensate is irrotational and thus can only 
acquire angular momentum through the nucleation of 
vortices [21]. For small rotation rates, the condensate's 
ground state has zero angular momentum. As Q increases 
and more vortices penetrate the cloud, the system under- 
goes a series of transitions from states with low to high 
angular momentum [22], eventually yielding incompress- 
ible QH-like states as J7 ^ 1. Fig. 1 shows the cor- 
responding 'yrast line' [23] for six particles. Currently, 
rotation rates in excess of f2 = 0.99 have been achieved 
in the laboratory [24] , though with current experimental 
conditions the filling factor v is still an order of magni- 
tude too high to reach the QH limit. 

The addition of V^xt has a profound effect on the 
RBG's energy spectrum since a narrow gaussian primar- 
ily increases the energy of levels with small angular mo- 
mentum L, because particles in these states arc situated 
closer to the origin. Consequently, including a repulsive 
Gaussian favors ground states with higher L, lowering the 
various critical rotation frequencies. As shown in Fig. 1, 
increasing 7 at constant cr, and fi, changes the ground 



4 - 



2 



3 



T ' 1 ' 1 ' 1 ^ 

2.0 

IL... 1-5 
:1.0 



0.0 



\ 

\ 

\ 



&H0.5- • ^. 
— [ ."— — 

"C^ 



L = 

- - - L = N 

Pfaffian 

Laughlin 



1 2 3 4 5 

p(0 




FIG. 2: Second order radial pair correlation function for the 
L — 0, L = N , Pfafhan, and Laughlin states of a six-particle 
RBG, with one particle fixed at the origin (p' = 0), in the 
absence of any external potential (7 = 0). The inset shows 
the corresponding densities. 

state from a Bose condensate with L = to the 'sin- 
gle vortex' state with L = N, eventually to the so-called 
Pfaffian state a.t L = N{N — 2)/2, and ultimately to the 
Laughlin state at L = N{N - 1) [22]. 

While higher angular momentum states may be fa- 
vored by the presence of Ve^t, and the resulting particle 
density might look similar as in its absence, it is not ob- 
vious that the strongly correlated nature of these states 
is preserved. To address this issue one needs to con- 
sider higher-order correlation functions. In its normal- 
ized form, the pair (also known as the density-density or 
second-order) correlation function is defined as: 



Vit(r)Vi^(r')V'(r')V'W 



^t(r')V,(r')) (V't(r)V)(r) 



(2) 



where V'^(r) and V;(r) are the bosonic creation and de- 
struction field operators, respectively, and the expecta- 
tion value (• • • ) is taken with respect to the ground state. 
The density- density correlation function can be inter- 
preted as the probability of a particle being at position 
r when another is at position r', and has been used to 
probe for pairing in a superfluid Fermi gas [25, 26] and 
for density-wave order in bosons confined within deep 
optical lattices [27-29]. 

As shown in Fig. 2, the Laughlin state has a particu- 
larly distinctive radial pair correlation function g2(p, p') 
(the (j) dependence is not accessible due to the cylindri- 
cal symmetry of our system), which would make its ex- 
perimental detection unambiguous. The Laughlin wave- 
function is an exact solution of the many-body Hamilto- 
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FIG. 3; Radial pair correlation for the five-particle Laugh- 
lin state, in the absence of any external potential (a). As 
the external potential amplitude 7 is increased to 7 = 0.25 
(b), 7 = 0.5 (c), and 7 = 1.0 (d), the likelyhood of overlap- 
ping particles increases suggesting that the ground state is no 
longer the Laughliri. Parameters correspond to g = 0.1, and 
a = 1.4. 




FIG. 4: Optimized magnitude (7) and width (a) of the exter- 
nal gaussian potential such that the RBG ground state is the 
Laughlin state for i}^ = 0.99 and g = 0.1. The lines are best 
fits, corres ponding to 7 = (0.039 ± 0.0 02) iV - (0.10 ± 0.01) 
and a = ^y{1.74 ± 0.06)Ar - (4.9 ± 0.3). 



nian (1) when 0=1, and is defined as 

*Laughii„ =Aff[{zi - z,)^e-^" \^^\'/<, (3) 

i<j 

where Zj = pj-e^'^J is the complex coordinate of particle 
j, and A/" is a normalization constant. The characteristic 
length £c = i/Ti/mWc is defined in terms of the cyclotron 
frequency ojc = 20, so that £c ~ £/V^ when O ^ cj. Be- 
cause of the Jastrow-type prefactor, the value of 52(1", r') 
is expected to vanish when two particles overlap (r=r'). 
In this respect the Laiighlin state is most different from 
the weakly correlated Bose-Einstein condensates at low 
L, for which 512 (r,r') is a constant approaching unity as 
N ^ 00 for any r and r'. 

Fig. 3 shows g2{p,p') with p, p' < 3£, for the ground 
state at O = 1 under the influence of different external 
potential intensities 7 and fixed a = 1.4. If t4xi had 
no influence on the Laughlin state, the whole diagonal 
would remain zero. Since this is evidently not the case, 
it is critical to choose the parameters 7 and a judiciously 
to best preserve the fundamental characteristics of the 
Laughlin state, particularly its pair correlations. 

For each N the parameters defining T4xi were varied in 
the range {0, 20} in their appropriate units, with a and 
7 in increments of 0.1 and 0.01, respectively. For each 
case, the location of the radial pair correlation maximum 
52(Pm,)Pm) along the diagonal was determined, and any 
parameters for which the maximum exceeded 0.01 were 
discarded. The values of the parameters 7 and a were 
minimized, subject to the condition that the RBG ground 
state is the Laughlin state for some chosen, fixed, value of 
the rotation frequency = 0/, < 1. Results are plotted 
in Fig. 4 for O^ = 0.99 and g = 0.1. These indicate that 



a linear increase in 7 with N and a ex a/Iv are sufficient 
to attain the Laughlin state at CIl = 0.99. Extrapolating 
from the small- iV data in Fig. 4, one obtains 7 = 39 ± 2 
and a = 42 ± 1 for N = 10^. Even for N = 10^, the 
values are not completely unreasonable: 7 w 400 and 
(J « 140. For a choice of a; = lOOHz, this corresponds to 
an amplitude of 7//1 = 40kHz, while optical lattices with 
depths of 50kHz are readily produced [16]. 

The qualitative behavior described above was found 
to be quite general. The numerical calculations were re- 
peated using two and three landau Levels, and the re- 
sults were identical to those obtained within the low- 
est Landau level approximation for interaction strengths 

< 5f < 1. The optimal values of the Gaussian pa- 
rameters depend on both g and O, but their functional 
dependence on N is unchanged. When g is increased, 
it is easier to reach the quantum Hall states [20], and 
therefore the values of cr and 7 are smaller. Similarly, 
should a lower transition rotation frequency O^ be de- 
sired keeping g constant, then the Gaussian needs to be 
higher and wider. For example, for O^ = 0.95 the best 
fits to the data yield 7 = (1.05 ± .05)iV - (2.3 ± 0.3) 
and a = ^(10.2 ± 0.5)iv - (23 ± 3) which yields extrap- 
olated values 7 = 1050±50 and a = 101±3 for N = 1000. 
Evidently, a very low transition frequency for the Laugh- 
lin state is not feasible with this technique. 

To further validate the pair correlation function as a 
figure of merit, one requires high overlap (^'Laughl^) be- 
tween the exact Laughlin state and the one obtained in 
the presence of the external potential. For all N con- 
sidered, the overlap integral exceeded 0.99. Plotting the 
overlap as a function of 1/N shows an asymptotic be- 
havior, with a limiting value for large N of 98.5% for the 

01 = 0.99 data shown in Fig. 4. Thus, for the optimized 
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values of 7 and cr, the state remains the bosonic fractional 
quantum Hall Laughlin state. 

A simple model accounts for the success of applied 
Gaussians in reducing the critical frequency for the 
Laughlin state: the combination of the external poten- 
tials and the centrifugal barrier should mimic a constant 
potential in the spatial region where the particle den- 
sity is appreciable. For large a, the Gaussian can be 
expanded in a Taylor series and the radial contribution 
to the single-particle potential V{pi) in Eq. (1) can be 
approximated as 

^^7+2 [1-^^ -^\^ + ^ + 0{Pi). (4) 

Enforcing that the term in the square brackets vanishes 
when Q = Ql, one obtains the relationship between the 
parameters of the Gaussian: 7 = (l — a^. To ensure 
that the Laughlin state experiences little of the quartic 
contribution in Eq. (4), one requires that the RMS width 
of the Laughlin state satisfy prms = \/ {p^) ^ cr for all 
N . The spread of \E'Laughiin can be estimated by not- 
ing that the Jastrow prefactor in the definition (3) can 
be written as a polynomial of order N{N — 1) in the 
variables pi with complex coefficients. Thus we obtain 
^Laughlin « A/"e'^(^-^)'»p^(^ -^)e-^p'. Stra ightforward 
integration yields prms = -y/iV — 1 + {1/N) w \/N for 

» 1. Indeed, direct evaluation using (3) for small A'^ 
gives exactly Prms = V^^- Enforcing the condition that 
Prms ^ cr, one obtains a ^ ^/N and 7 ^ A/', in agreement 
with the results shown in Fig. 4. 

The best fits to the numerical data are consistent with 
these scaling relations. For g — 0.1 and = 0.99, 
one obtains a = \/1.74:N — 4.9, yielding an optimal 
value for the amplitude of 7 = (1.747V - 4.9)(1 - 
from which the correct slope and y-intercept values are 
found. Intriguingly, for fi/, = 0.95, one rather obtains 
a = \/lQ.2N — 23, so that the Gaussian width must be 
considerably wider than the mean spread of the Laughlin 
density as the desired rotation frequency is lowered. It is 
important to note that these choices for the Gaussian pa- 
rameters facilitate access to not only the Laughlin state, 
but to all ground states: the effective trap frequency is 
reduced from \/\ — ^l? to -^O^ — Q?, so all transition 
frequencies are accordingly reduced. 

In conclusion, we have demonstrated that it is possi- 
ble to attain the Bose-Laughlin quantum Hall state, and 
thus any other quantum Hall state, by applying an exter- 
nal Gaussian potential to a harmonically trapped rotat- 
ing Bose gas. With experimentally reasonable choices of 
Gaussian amplitude and width, number of particles, and 
effective interaction strength, the transition to a Laughlin 
state can occur for the experimentally accessible rotation 
rate of 0.99a; or lower. These results suggest that quan- 
tum Hall physics in neutral, rotating Bose systems are 
currently within experimental reach. 
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